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Abstract — "Extended Clifford algebras" are introduced as 
a means to obtain low ML decoding complexity space-time 
block codes. Using left regular matrix representations of two 
specific classes of extended Clifford algebras, two systematic 
algebraic constructions of full diversity Distributed Space-Time 
Codes (DSTCs) are provided for any power of two number of 
relays. The left regular matrix representation has been shown 
to naturally result in space-time codes meeting the additional 
constraints required for DSTCs. The DSTCs so constructed 
have the salient feature of reduced Maximum Likelihood (ML) 
decoding complexity. In particular, the ML decoding of these 
codes can be performed by applying the lattice decoder algorithm 
on a lattice of four times lesser dimension than what is required 
in general. Moreover these codes have a uniform distribution of 
power among the relays and in time, thus leading to a low Peak 
to Average Power Ratio at the relays. 

I. Introduction 

Coding for wireless relay networks has received a lot of 
attraction recently with the advent of cooperative diversity 
techniques. In this paper, we are interested in constructing 
Distributed Space-Time Codes (DSTCs) for the Amplify and 
Forward (AF) based cooperative diversity protocol proposed 
by Jing and Hassibi [1]. The Jing and Hassibi protocol is a 
two phase based AF protocol. In the first phase, the source 
broadcasts a vector to all the R relays which contains the 
information that the source intends to communicate to the 
destination. In the second phase, each relay transmits a vector 
obtained by linear processing of the received vector and its 
conjugate to the destination. To the destination, this would 
appear as if each relay transmitted a column of a linear space- 
time code thus leading to the concept of DSTCs. We refer the 
readers to [1] for a detailed introduction to DSTCs. Consider a 
RxR linear design S(X) in 2K real variables x\, x%, . . . , x%k 
as follows 

IK 

S{X) = Y d X i C i 

z=l 

where, the complex matrix Cj <E c«x-R j s called the 'weight 
matrix' corresponding to the real variable Xi and R is the 
number of relays. Let X = [ x\ X2 ■ ■ ■ x-ik ] ■ From 
the 2K real variables we can form K complex variables 
zi,Z2,--- ,Zk by pairing two real variables at a time. Let 
s = [ z\ Z2 ... zk ] • Then the linear design S can 



also be expressed as 

S{X) = [ Ais + B lS * A 2 s + B 2 s* ... A R s + B R s*] 

(1) 

where, Ai, Bi 6 C RxR , i = 1, . . . , R are complex matrices. 
We call the matrices Ai,Bi as 'relay matrices'. In [1], it 
has been shown that linear designs satisfying the following 
conditions 

1 ) For any i — 1 , . . . , R, either Ai = or B, = 
(Conjugate-Linearity Property) and K = R 

2) All the nonzero relay matrices are unitary matrices 

are applicable as DSTCs where s = [ Z\ Z2 ... zr ] T 
will be the vector transmitted by the source in the first phase 
and the matrices Ai , Bi will be used at the ith relay to perform 
linear processing of the received vector and its conjugate. 

All the previous works on DSTC construction except [2], 
[3], [6] and [7] do not address the important problem of 
designing DSTCs with low Maximum Likelihood (ML) de- 
coding complexity. This problem gains significant importance 
especially if the number of relays in the network is large. 
Suppose we partition the 2K real variables and their corre- 
sponding weight matrices into ^-groups Lk,k = 1, . . . ,g, the 
fc-th group containing 2K/ g real variables. Then without loss 
of generality we can consider the natural simplest partition, in 
terms of which S(X) can be written as, 

IkK 

9 3 

s(x) = y. s *( x ^ where > s k(Xk)= J2 XiCi - 



If the received matrix and the channel matrix are denoted by Y 
and H respectively, then a ML decoder in general minimizes 
the metric || Y ~ S(X)H || 2 . However, if 

Cf Cj + CfC t = 0, V^ e L p , Vj ei„p/?,l< M <s 

(2) 

then this is equivalent [6] to minimizing 

|| Y-S k (X k )H || 2 (3) 

for each 1 < k < g individually. Note that (O can be computed 
by applying the lattice decoder algorithm on a lattice of g 
times lesser dimension. We then say that the DSTC is g-group 
ML decodable. In [2], the authors constructed 2-group ML 



decodable codes using division algebras. In [3], DSTCs for 
two and four relays based on the Alamouti design and quasi- 
orthogonal design were proposed. Following the works of [2], 
in [6], a class of 4-group ML decodable codes called Precoded 
Co-ordinate Interleaved Orthogonal Designs (PCIODs) were 
constructed for arbitrary number of relays. However, PCIODs 
have a drawback that the power distribution among the relays 
is not uniform across time slots thus leading to a large Peak to 
Average Power Ratio (PAPR). Moreover the relay matrices of 
PCIODs are not unitary which in turn forces the destination 
to perform additional processing [6] to make the covariance 
matrix of the resultant noise vector at the destination a scaled 
identity matrix. Recently in [7], this problem was resolved by 
an alternative iterative construction for number of relays R 
being a power of two. However, both the constructions of [6], 
[7] were not obtained from a systematic algebraic procedure 
targeting the requirements for low decoding complexity. The 
main contributions of this paper are as follows. 

• A generalization of Clifford algebras, which we call "Ex- 
tended Clifford Algebras" is introduced as an algebraic 
framework to handle the problem of constructing STBCs 
with low ML decoding complexity. To the knowledge of 
the authors, this is the first known systematic algebraic 
procedure to solve this problem. This algebraic frame- 
work simplifies the problem to finding appropriate matrix 
representations of extended Clifford algebras. 

• Using left regular representation of "Extended Clifford 
Algebras", two different fully diverse algebraic DSTC 
constructions are provided for power of two number 
of relays. Left regular representation has been shown 
to naturally result in space-time codes meeting the ad- 
ditional requirements of DSTCs. Moreover, one of the 
constructions provides an algebraic explanation for the 
recently proposed DSTC construction in [7]. 

• ML decoding of these algebraic DSTCs can be performed 
by applying the well known lattice decoder algorithm 
on a lattice of four times lesser dimension than what is 
required in general. 

• Furthermore, the proposed DSTCs have lesser PAPR than 
the DSTCs of similar ML decoding complexity reported 
in [6]. 

The rest of the paper is organized as follows: In Section HIl 
we briefly recollect a set of known sufficient conditions for low 
ML decoding complexity designs. Extended Clifford algebras 
are introduced in Section [Hi] and algebraic construction of 4- 
group ML decodable DSTCs is presented in a general setting. 
The two special classes of codes from extended Clifford 
algebras are presented in detail in Sections |IV] and [V] In 
Section [VI] it is shown that all the requirements for DSTCs 
are satisfied by the proposed codes. 

Notation: For a complex matrix A, Aj denotes the real matrix 
obtained by taking the real parts of all the entries of A and 
Aq denotes the real matrix obtained by taking the imaginary 
parts of all the entries of A. If A is an algebra over a field F 
then EndF(A) denotes the set of all maps from A to A that 



are F— linear. 

II. Sufficient conditions for low ML decoding 

COMPLEXITY DESIGNS 

Recently in [4], sufficient conditions for designing g-group 
ML decodable STBCs have been reported. Since our con- 
structions rely upon these sufficient conditions and because 
we restrict ourselves to the g = 4 case, we briefly introduce 
these sufficient conditions [4] for g — 4 before proceeding 
further. Essentially we would like to be able to partition the 
set of weight matrices of a linear design into 4 groups in such 
a way that the condition in (fJJ) is satisfied. Let us first list 
down the K weight matrices in the form of an array as shown 
below. 







Ck _i_ i C 3K j_i 

, -1-1 4 -1-1 


c 2 


Ck +2 




Ck 

4 


Ck 

2 


C 3K Ck 

4 



The partitioning is as follows: All the weight matrices in one 
column belong to one group. To simplify the construction, 
we shall consider all the weight matrices to be unitary and 
furthermore set C\ — I. Then it has been shown in [4] that 
it is sufficient to design the matrices in the first row and the 
first column such that they satisfy the following conditions. 

1) All the matrices in the first row except C\ = I should 
square to — / and should pair- wise anti-commute among 
themselves. 

2) The matrices in the first column should square to / and 
should commute with all the matrices in the first row 
and the first column. 

Once such a set of matrices is obtained, the matrix in the i-th 
row and j-th column can be filled up by multiplying Ci and 
C (j-i)K , -i . It can be easily verified that such a set of weight 
matrices will satisfy the conditions in (O for g = 4. 

III. Algebraic construction 4-group ML 

DECODABLE DSTCS 

An algebra is simply a ring as well as a vector space 
with the addition operation being compatible to both the ring 
and the vector space structures. In this section, we introduce 
the algebraic framework of "Extended Clifford Algebras" to 
handle the problem of constructing g-group ML decodable 
codes satisfying the sufficient conditions discussed in the pre- 
vious section [7]. Using left regular representation of extended 
Clifford algebras we then obtain two constructions of 4-group 
ML decodable DSTCs. 

Our methodology to construct the matrices in the first row 
and first column (as discussed in previous section) would be 
to fabricate an algebra in such a way that it contains elements 
satisfying the algebraic relations we need. Once we construct 
the algebra, we then obtain the required linear design by 
taking an appropriate matrix representation of the constructed 
algebra. 



Definition 1: [5] The Clifford algebra, denoted by Cliff n 
is the algebra over the real field K generated by n ob- 
jects 7^, k — 1, ...,n which are anti-commuting (-fkjj = 

~7j7fc, Vfc 7^ j) and squaring to —1 (7! = — 1 Vfc = 
l,...,n). 

A natural basis for Clif ' f n seen as a vector space over M is 
^n= {l}U{7i|i = !,•••,»} m 

u;; l=2 {nr:i7 fcl ii<fc,<fc, + i<«} 

The number of basis elements is |^„| = 2". 

Notice that the defining algebraic relations of the generators 
of a Clifford algebra resemble the algebraic relations which 
the matrices in the first row should satisfy. Hence we can 
obtain the matrices in the first row by taking unitary matrix 
representations of the generators of a Clifford algebra. To ob- 
tain the matrices in the first column, we use a similar strategy. 
We introduce few new symbols in the Clifford algebra and 
define them to square to 1 and commute with the generators of 
the Clifford algebra and also commute among themselves. In 
other words, after introducing new symbols, multiplication in 
the algebra is appropriately defined in order to create a bigger 
algebra which contains the Clifford algebra as a sub-algebra. 
Hence by taking the unitary matrix representation of these 
specific elements of the algebra, we get the weight matrices 
of the required linear STBC. 

Definition 2: Let L = 2 a ,a e N. An Extended Clifford 
algebra denoted by A^ is the associative algebra over R gen- 
erated by n + fl objects k = 1, . . . , n and 5^ i = 1, . . . , a 
which satisfy the following relations: 

. 7 2 = -l, Vfc = l,...,n 

• 7fe7j = -7j'7fci V k ^ j 
. 5£ = 1, Vfc = 1,. . . ,a 

• S k Sj = S 3 S k , V 1 < k, j < a 

• <5fc7j = 7j^fc, Vl<Ha,l<j<Ti 

From the above definition, it is clear that Cliff n is a sub- 
algebra of A„. Let 38 n be the natural K basis for this sub- 
algebra Cliff n . Then a natural R basis for A^ is given by 

S%= 3§ n U {S8 n kV = 1, ■ ■ • , a} 

IC=2 ®n {Iltl ffc.|l < *i < h+1 < a} . (D) 

Thus the dimension of A,^ seen as a vector space over R is 

2'M+a 

Example 1: Let us take n = 2, a = 1. Hence L = 2. Then 

Aj = {ax + 7ia 2 + <5ia 3 + (5i7ia 4 |ai, a 2 , a 3 , a 4 e K} . 

Addition in the algebra is defined to be component wise and 
multiplication is completely described by defining the multi- 
plication between any two basis elements. The multiplication 
table can be easily generated using the defining algebraic 
relations of the generators and is given as follows . 





1 


71 


Si 


£171 


1 


1 


7i 


Si 


£171 


71 


71 


-1 


(5i7i 


-Si 


Si 


Si 


£171 


1 


7i 


£171 


£171 


-Si 


7i 


-1 



One can check from the multiplication table that the multipli- 
cation is indeed associative. Note that A\ can also be viewed 
as a vector space over C by thinking of the symbol 71 as the 
complex number i = \/—\. Then, we have 

A 2 = {zi + 5iz 2 \z u z 2 e C} 

where, zi = ai + 7102 and z 2 = a 3 + 7104 

Since we are interested in 4-group decodable DSTCs, we 
need 4 matrices (including identity matrix) in the first row. One 
way to obtain such matrices is to take the matrix representation 
of A3 for L = 2 a ,a € N. The matrix representation of the 
symbols l,7i,7 2 ,73 respectively can be used to fill up the 
first row. Interestingly, there is yet another way of obtaining 
such matrices. Let us look at A 2 for L = 2 a ,a 6 N. The 
symbols 71 and 72 square to —1 and anticommute. However 
note that 

(7271) 2 = —i; (7271)71 = -71(7271); (7271)72 = -72(7271)- 

Thus the symbol 7 2 7i also squares to —1 and anticommutes 
with the symbols 71 and 72. Thus we can fill up the first row 
with the matrix representations of the symbols 1,71,72,7271 
respectively. Thus we get two classes of 4-group ML decod- 
able STBCs, one from A3 and the other from A^ . 

A. Matrix Representation 

There are several ways to obtain a matrix representation of 
an algebra. However we need to take an appropriate matrix 
representation such that the following conditions are satisfied. 

1) The symbols 1, 71, 72, 7„, 5 k ,k = l,...,a, 
Um=2llj=i£fcil 1 < h < fci+i < a should be repre- 
sented by unitary matrices. 

2) The resulting linear design should have the Conjugate- 
Linearity property. 

3) All the relay matrices should be unitary. 

Such matrices are naturally provided by the left regular 
representation of the associative algebra A^. Left regular rep- 
resentation is an easy way to obtain the matrix representation 
for any finite dimensional associative algebra [8]. The first 
requirement of unitary matrix representation is met because the 
natural basis elements of A^ together with their negatives form 
a finite group under multiplication. This fact in conjunction 
with the properties of left regular representation guarantee a 
unitary matrix representation for the required symbols. We 
shall prove the other properties in Section [VT] after illustrating 
the construction procedure for both the codes from A2 as well 
as those from A3 . 

IV. Codes from A^ 

We first view A 2 as a vector space over C by thinking of 
71 as the complex number i = \/—l. A natural C basis for 
A^ is given by 

B£= {l,72}U{{l, 7 2}£ i |i=l,...,a} 

U™=2 {1, 72} {UZi S ki |1 < h < k l+1 < a} 



The dimension of Aj seen as a vector space over C is 2™+ a ~ 1 . 

We have a natural embedding of into Endc(A2 ) given 
by left multiplication [8] as shown below: 

: Aj h-> End c (A2 ); <f>(x) = L x : y i-> xy. 

Since the map L x is C -linear, we can write down a matrix 
representation of L x with respect to the natural C basis B^. 
Thus we obtain a design satisfying the requirements of (O for 

9 = 4- 

Example 2: Let us begin with Nt — 2 transmit antennas. 
Let n = 2. Then equating n + a — 1 = 1, we get a = 
and hence L = 1. But the algebra A3, is same as Cliff (2) 
which is nothing but the Hamiltonian Quaternions EL It is 
well known [8] that the left regular matrix representation of 
H yields the popular Alamouti design. Thus we see that our 
algebraic code construction which was driven by the need for 
low ML decoding complexity naturally leads to the Alamouti 
design. 

Example 3: Suppose we want a design for Nt = 8 = 2 3 
transmit antennas. Let n — 2. Then we need n + a — 1 = 3. 
Thus a — 2 and L — 4. A general element of the algebra Af 
looks like 

X = 21 + SlZ2 + #223 + 8l52Z4 + 72Z5 + #17226 + #27227 + 5l8272Z8 

where, Zi € C, Vi = 1, . . . , 8. The image of the basis £>| under 
the map is shown in © at the top of the next page. Thus, 
we have 



Zl 


22 


23 


24 


-25 


-26 


-2? 


-28 


22 


Zl 


24 


23 


-26 


-25 


-2g 


-2? 


23 


24 


21 


22 


-2? 


-zl 


-25 


-26 


Z4 


23 


22 


Zl 


-2l 


-2? 


-26 


-25 


25 


26 


27 


28 


z* 


22 


23 


24 


26 


25 


28 


27 


22 


2* 


zl 


4 


27 


28 


25 


26 


23 


24 


Zl 


22 


28 


Zl 


26 


25 


2 4 ' 


2j 


22 


2l - 



Also, we have 

X = 21/ + 7l2lQ + 5l2 2 J + 5l7l2 2 Q 

+52231 + <527l 23Q + 8182Z41 + 8lS 2 'yiZ4Q 
+7225/ + 7271 25Q + #17226/ + <5l727l26Q 
+#27227/ + #27271 27Q + #1#27228/ + #l#2727l28Q 

leading to 

L x = 0(1)21/0(1) + <f>(ji)ziQ + 4>{5 1 )z 2 i 

+0(#l7l)2 2 Q + (p(5 2 )z 3I + 0(#27l) 2 3Q 
+0(<5l#2)24/ + 0(#l#27l)24Q +0(72)25/ 
+0(727l)25Q + 0(#l72)26/ + 0(#l727l)26Q 
+0(#272)27/ + 0(#2727l)27Q + 0(#1#272)28/ 
+0(#l#2727l) 2 8Q 

which explicitly gives the design L x in terms of its weight 
matrices. Expressing the elements of the algebra, the real 
variables of the resulting design and their corresponding 
weight matrices in the form of a tabular column as discussed 
in Section UT1 we get 



1 


71 


72 


7271 




0(7i) 


0(72) 


0(727l) 


Zll 


ZlQ 


25/ 


25Q 


Si 


Slfl 


#172 


#17271 


HSi) 


0(5i)0( 7 i) 


0(#l)0(72) 


0(#i)0(727i) 


22/ 


22Q 


26/ 


26Q 


<5 2 


<527l 


#272 


#27271 


4>{S 2 ) 


^(#2)0(71) 


^(#2)^(72) 


0(#2)0(727l) 


23/ 


23Q 


27/ 


27Q 


&1&2 


(5ii527i 


#1#272 


#1 #27271 


HS1S2) 


(f>(5i #2)0(71) 


0(#1#2)0(72) 


0(#l#2)0(727l) 


Z4J 


2 4 Q 


28/ 


28Q 



Lrom the table above, it is clear how the weight matrices 
and real variables can be partitioned into four groups. 



In general for Nt — 2 transmit antennas we take the 
left regular representation of A2 to obtain a 4-group ML 
decodable linear design satisfying (ffji for g = 4. These codes 
were first obtained using a non-algebraic iterative construction 
procedure in [7]. The algebraic framework presented here 
provides an interesting algebraic explanation for the codes in 
[7]. 

V. Codes from A3 

We use a slightly different approach to obtain codes from 
A3. Let us first consider the algebra, A3 which is nothing but 
Cliff 3. A general element of Cliffs looks like 

x = ai+7ia2+72«3+73«4+7i 7205+7273^6+71 73«7+7i7273a8 

for some a, € R, i = 1, ... ,8. The element 717273 satisfies 
the following properties. 

(717273) 2 = 1; (71X717273) = (717273X71); 
(72) (717273) = (717273) (72); (73) (717273) = (717273) (73)- 

Thus the element 717273 squares to 1 and commutes with 
all the generators of Cliffs. Hence the matrix representation 
of the element 717273 can be used as a candidate to fill up 
the first column. Since we have now filled up two matrices 
(including the identity matrix) in the first column, it should be 
possible to get a 2-real symbol decodable code using matrix 
representation of Cliffs. From Section HU we know that the 
remaining weight matrices should be obtained as a product 
of matrices in the first row and those in the first column. We 
have, 

(7l)(7l7273) = -7273; (72)(7l7273) = 7i73; (73)(7l7273) = -7i72- 

It so turns out that the elements 

{I,7i,72,73,-7i72,-7273,7i73,7i7273} also form a 
basis for Cliffs. Thus a general element of Cliffs can be 
expressed as 

x = ai + 7i«2 + 7203 + 7304 

+ (-7l72)a5 + (-7273)06 + (7l73)a7 + 71727308 

for some ai € R, i = 1, . . . , 8. By thinking of the element 71 
as the complex number i = we can view Cliffs as a 

vector space over C. To be precise, 

x = (ai +71^2) +72(03 +7105) 

+73(04 - 7107) + 7273(-a 6 + 71 a 8) 
= zi + 7222 + 7323 + 727324 



</>(!) 

0(72) 

0(5172) 
0(^72) 

<K#i#272) 



zi + 5 1 z 2 + S 2 z 3 + S 1 6 2 z 4 + 72^5 + #172^6 + #27227 + #i#272£8 
Sizi + z 2 + Si5 2 z3 + 5 2 Zi + #17225 + 72^6 + Si5 2 j 2 z 7 + S 2 ^ 2 z 8 
S 2 zi + 5\5 2 z 2 + 23 + SiZi + 8 2 ~/ 2 z 5 + 5i5 2 ~f 2 z e + 72 2:7 + #17223 
S 2 S 2 zi + 5 2 z 2 + 8iz 3 + z 4 + 5i6 2 ~/ 2 z 5 + 5 2 ~f 2 z e + #17227 + -y 2 z 8 
(zi + 5iz 2 + 5 2 z 3 + #i# 2 z 4 + 72Z5 + #17226 + £272^7 + 5iS 2 ^f 2 zs) 72 
J2Z* + <5i72z| + #27223 + <5i<5 2 7 2 24 - z% ~ hzt - #227 - #i# 2 2| 



(6) 



#172^ + 724 + 5ifal2zl + #272^4 - S 1 Z 5 



5\5 2 z 7 



5 2 zl 



faliZ* + #i#2724 + 72Z3 + #1722:4 - #225 - <5 2 26 - z? - SiZg 

+ &2l2Z* 2 + #1724 + 7 2 2| - #1#225 - #22| _ &l z 7 ~ z 8 



where, Zi € C, i = 1, . . . ,4 and are given by 

zi = (ai + 7102); Z2 = (03 + 7105); 
23 = (04 - 7i fl 7); Z4, = (-a 6 + 7ia 8 ). 

Now using left regular representation as in the case of codes 
from A2 , we obtain the following design 
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In general, for R = 2 A relays we take the left regular 
representation of A3 

Example 4: Suppose we want a design for R = 8 = 
2 3 relays. Hence we have A = 3. Using the left regular 
representations of the algebra A§, we get the following linear 
design 



Zl 


-zl 


-zi 


-Zi 


Z5 


-z% 


-z? 


-zs 


Zl 


zl 


-z 4 


Z3 


Z6 


zt 


-zi 


Z7 


Z3 


zl 


z* 


— Z2 


ZT 


zt 


zi 


-ze 


Zi 


-zi 


zl 


Zl 


Z8 


-Zl 


zi 


zs 


Z5 


-zl 


-Z7 


-zs 


Zl 


-zl 


-zi 


-Zi 


Z6 


zi 
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-zl 


zl 
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The corresponding 4 groups of real variables are 

{Zll, Z4Q1 Z5I, Z S q}, {ZIQ, Z41, 2 5 Q, Zgl}, {z 2 I, 2 3 Q, Z 6I , Z7q} 

and {z 3I ,z 2 q,z 7I ,z 6Q }. 

VI. Conjugate-linearity, unitary relay matrices 

AND FULL-DIVERSITY 

Note that both the classes of codes from Aj and A3 have the 
property that any column of the design has only the variables 
or their conjugates. This is by virtue of the properties of left 
regular representation. While taking the left regular matrix 
representation, recall that we viewed the algebra as a vector 
space over C by thinking of the element 71 as the analogue of 
the complex number i = ^— 1, Any column of the design 
was then obtained as the image of a few elements of the 
natural basis of the algebra under the map L x . All the elements 
of the natural basis of A^ have the property that they either 
commute with 71 or anticommute with 71. When we found the 
image of a basis element say y, recall that we moved y past 
a complex number Zj. If y commutes with 71, then it leaves 



the complex number intact. If y anticommutes with 71, then 
it inflicts conjugation while moving past the complex number. 
This fact can be clearly observed in ©. 

Moreover, it can be easily observed that all the relay 
matrices of the resulting designs are unitary. This is because 
the number of complex variables in the design is equal to 
the size of the matrix and by virtue of the left regular 
representation any complex variable appears only once in any 
column. Further, the positions in which they appear in different 
columns is different. 

Full diversity can be obtained for all the constructed codes 
by choosing an appropriate rotated lattice constellation. 
This has been proved in [4] more generally for all the 
codes constructed using the sufficient conditions for low ML 
decoding complexity discussed in Section Ull 
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